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ABSTRACT 

We propose a new mechanism for the delivery of gas to the heart of galactic nuclei. We show 
that warm halo clouds must periodically impact galactic centers and potentially deliver a large (~ 
10 4 ~ 6 Af Q ) mass of gas to the galactic nucleus in a singular event. The impact of an accreting warm 
halo cloud originating far in the galactic halo can, depending on mixing, produce a nuclear starburst 
of low metallicity stars as well as low luminosity accretion onto the central black hole. Based on 
multiphase cooling around a ACDM distribution of halos we calculate the nuclear impact rate, the 
mass captured by the central black hole and the fraction of active nuclei for impacting cloud masses 
in the range 10 4 — 10 6 M Q . If there is moderate braking during cloud infall, our model predicts an 
average fraction of low luminosity active nuclei consistent with observations. 



1. INTRODUCTION 

Activity in galactic nuclei is fuelled by a reser- 
voir of low angular momentum gas, but it is unclear 
how such reservoirs build up. Fast outflows from 
OB stars, su pernovae or AGN activ i ty can clear out 
nuclear ISM (ICiotti fc Ostrikerl Il997t iSchawinski et al.l 
120071 iMcKernan et al.l l2007t iSchartmann et al.l I2Q09TT 
Gas in the nucl eus can also be quickly consumed 
by star formation ( Navaksh in. Cuadra fe Springell 120071) 
or driven outwards by p ositive gravitational torques 
(jGarcia-Burillo et al.l I2005T) . If the gas in the reservoir 
originates outside the nucleus, multiple mechanisms op- 
erating on different distan ce- and time -scales, such as 
bars feeding nuclear rin gs (|Krohkll2001l ) radiation drag 
dKawakatu et al.l 12003?) or interactions (jVittorini et al.l 
I2005D are required. From observations, around 1/3 of 
all galactic nuclei in the l ocal Univ erse exhibit low lu- 
minosity nuclear activity (Ho 2008), so the mechanism 
promoting nuclear gas build up is likely to be simple and 
on-going. 

In this letter we propose a new mechanism for delivery 
of gas to the galactic nucleus. We show that the impact 
of a warm halo cloud (WHC), containing - 1O 4 " 6 M , 
on the central regions of a galaxy will fuel nuclear activ- 
ity. We develop a simplistic model of this phenomenon to 
demonstrate its likely importance. While there remains 
considerable uncertainty in WHC parameters, we show 
that for plausible input parameters, a direct hit by a sin- 
gle WHC on the center of a galaxy will supply fuel for star 
formation and radiatively inefficient accretion onto the 
central black hole. Since WHC bombardment of galaxies 
must occur, some fraction of the low luminosity activity 
observed in galactic nuclei must be due to WHC impacts 
on galactic nuclei. 

2. WARM HALO CLOUDS 
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The assembly of gaseous halos around galaxies nat- 
urally produces a multiphase medium with warm halo 
cloud s (WHCs) embedded i n a low density hot gas 
halo (IMaller fc Bullock! l200l . WHCs could be the lo- 
cal analogue of t he high redshift Lyman limit systems 
(jBlitz et al.l 119991) . In the halo of our own Galaxy, nu- 
merous warm clouds are observed with vel ocities that 
are in consistent with Galactic rotation (e.g. Pe ek et all 
2007). Around several other galaxi es, HI clouds of 
\Q b ~ s M G) have been detec ted (e.g. iChvnoweth et al.l 
120081: iBoomsma eFaXl I2008D . Extended HI structures 
(up to ~ 10 10 Mq) are observed out to ~ lOOkpc around 
early- type galaxies (lOosterloo et al.l l2007). So, a popula- 
tion of WHCs containing ~ 1O 4_6 M0 per cloud, may be 
comm on around most galaxies (see e.g. IMaller fc Bullock! 
2004). Around our Galaxy, two basic models can ac- 
count for the properties of observed high velocity clouds 
(HVCs): an ac c retion model based on WHCs (e.g. 
IBlitz et al.l [19991: IMaller fe Bullockl l200l or a Galac- 
tic fountain model (Sha piro fc Field! I1976D . Accreting 
WHCs should dominate the mass of clouds in the Galac- 
tic halo so in this work we shall concentrate on the effects 
of an accreting WHC impact on galactic nuclei. Much 
of the following discussion also applies to galactic foun- 
tain clouds, although that population will have higher 
metallicity, less mass, smaller radius and lower velocity 
on average, and will be less numerous. 

3. WHC IMPACT WITH GALACTIC CENTER 

Large uncertainties exist concern ing cloud trajecto - 
ries around our own Galaxy (e.g. iThom et al.l [20081. 
While we might naively expect radial trajectories for 
WHCs, cloud trajectories may become randomized close 
to a galax y. Clouds close to the disk become tid ally 
disrupted (Blitz et al.l Il999t IMaller fc Bullockl I2004D. or 
deflec ted by magnetic fields (jKwak. Shelton fc Ralevl 
2009), or disrupted by the Kelv in-Helmholz instability 
(KHI) (jHeitsch fc Putmanl[2009T l and the fragments dis- 
persed about the disk, whereupon the previous cloud tra- 
jectory may be irrelevant. Since we have no way of pre- 
dicting actual cloud trajectories, in the discussion below 
we simply assume that clouds impact the disk randomly. 

We start by assuming a initial population of N c i clouds 
on random trajectories raining down on a galaxy of ra- 



2 



McKernan,Maller & Ford 



dius Rgai- As long as the WHC radius is larger than the 
nuclear region under consideration, the rate of impact 
(dNd/dt) of infalling WHC of radii r c i on the galactic 
center is 

dNd/dt = 0.25(r cl /R gal ) 2 / T (1) 

where r is the typical cloud infall time and the cloud 
material is assumed to arrive within 0.5r c ; of the galac- 
tic center. We calculate the cloud impact rate for mul- 
tiphase cooling aroun d a ACDM distribution of halos 
based on the model of IMaller fe Bullockl (|2004 ). In this 
model the total mass in WHCs is based on the mean 
free path of the clouds and the cloud properties we 
use are the average of some distribution in the halo. 
First we calculate the average nuclear impact rate, us- 
ing R ga i = 0.15f? eoo ;, based on con servation of angular 
momentum (Fall & Efst athioul|1980D. where R coo[ is the 
cooling radius (|White fc Fre nk 1991). We simplify our 
calculation by assuming that the WHC does not frag- 
ment and that infall time is r = R CO ol /V max where V max 
is the maximum circular velocity in the halo. The results 
are shown in Fig. [IJa), where we plot impacts/Gyr ver- 
sus Vmax for impacting cloud masses of 10 4 , 10 5 , 10 6 M Q . 
Evidently the impact rate is relatively flat with V max (or 
equivalently black hole mass). Note that a nuclear im- 
pact rate of ~ 5/Gyr from Fig. [IJa) corresponds to a 
galactic impact rate of ~ 200 / Gyr or a mass inflow rate 
of ~ lM /yr if the typical impactor mass is ~ 10 6 Af Q . 
This is approximately the low metallicity inflow rate re- 
quired to expla in the so-called 'G-dwarf problem' (e.g. 
Edmunds 1990, and references therein) . 

4. DELIVERING MATERIAL TO THE GALACTIC 
NUCLEUS 

Although obser vational constraints o f cloud impacts 
are not strong (|Lockman et al.l 120081 ). we can con- 
sider the consequences of WHC impact on a galac- 
tic nucleus, guided by simulations of HVC impacts 
with the disk (see e.g. iKwak. Shelton fc Ralevl 120091 : 
IHeitsch fc Putmanl 120091 and references therein). The 
same basic sources of cloud fragmentation in the halo and 
disk will apply to a WHC falling into a galactic bulge. 
If galactic nuclei have r elatively strong mag netic field 
strengths in general (e.g. iMetzger et al .1 119 96;). the mag- 
netic fields will act as a strong brake on infalling WHCs 
and the cloud may fragment. 

Instabilities in the infalling cloud have growth times 
of the order of the timescale on which shocks cross the 
cloud. This cloud crushing t imescale, t cc , is given by 
(IKlein. McKee fe Colellalll99l 
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n c i 
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where n c i and the densities of the WHC and the 

surrounding medium respectively and v c i is the relative 
velocity of the WHC and the surrounding medium. If the 
WHC is accreting from a very large distance, t cc can be 
less than the infall time (~ Gyrs) and the cloud can frag- 
ment in the halo. Indeed many of the HVCs presently 
observed may be fragments of initially much larger ac- 
creting WHCs. Galac t ic haloes are hot and d iffuse (e.g. 
iStrickland et alj 12004 iMcKernan et al.l 12004 and refer- 
ences therein) and the disruption of a WHC in the halo 



could create an infalling stringy association of clouds like 
the HVC A and C complexes. If the trajectory of such 
a stringy cloud complex intercepted the galactic nucleus 
then most of the WHC mass could be delivered to the 
region around the central supermassive black hole. 

Sufficiently dense and fast clouds falling through the 
halo and then the bulge ISM will first experience de- 
celeration and compression (pancaking), followed by an 
expansion of the shocked cloud downstream, then lat- 
eral expansion and finally cloud destruction when in- 
stabilities and differential forces fragment the cloud. 
This entire process takes plac e over a few t cc (see e.g. 



Klein. McKee fe Colellal Il994 IKwak. Shelton fe Ralevl 
2009f h Compression wi ll increase the WHC density 



as the cloud pancakes (|Kwak. Shelton fc Ralevl 120091 : 
IHeitsch fc Putmanl 120091 ). Small, fast, dense fragments 
of cloud ca n also sweep up large amoun ts of gas and 
be shocked (jWakker fc van Woerden|[l997D . The transit 
time for a WHC through a bulge of radius Rb to the 
galactic center is ~ Rb/rjv c i where 77 = [0, 1] is a coeffi- 
cient incorporating slowdown due to drag and magnetic 
fields. For a galaxy the size of our own(i?b ~ 3kpc), a 
cloud with v c i ~ 100 km s _1 will take ~ 30/?yMyr to 
cross the bulge and reach the center. 

But will the cloud survive this long? Outside of the in- 
ner ~0.5kpc, the bulge IS M in our own galax y has very 
low density (< !Q- 3 cm- 3 ) Prang et al.l (|2009f) . At these 
densities, the cloud crushing timescale is > lOOMyrs, far 
larger than the transit time through this part of the 
bulge. Within 0.5kpc of the Galactic center th e bulg e 
ISM density increases to ~ 0.01cm" 3 (|Tang et al.ll2009n . 
Here the crushing timescale is ~ 30/?7Myr, short enough 
that the WHC will fragment in this region. So, WHCs 
could survive until they reach the central regions, where- 
upon fragmentation is likely. 

5. GRAVITATIONAL CAPTURE BY THE CENTRAL 
BLACK HOLE 

Evidently a WHC cloud impact could deliver a large 
quantity of moderate density gas to a galactic nucleus in 
a single cloud infall. But what happens once fragments 
of a cloud of mass M c i ~ 10 4 ~ 6 M Q arrives in the central 
few tens of pc of a galaxy? Low angular momentum 
fragments of the WHC will be captured by the central 
black hole out to a radius (R eS c) where their velocity is 
below the escape speed. Thus the mass captured is 



R 2 

M capt = M c i-f^ 

T cl 



M, 



(2GM_ 



BH) 



(3) 



We calculate the average M capt for clouds of mass 
10 4_6 M Q impacting galactic nuclei. We estimate 
black hole mass usi ng the Mbh ~ a relation from 
iGultekin et al.l ()2009f ) and we assume that the bulge ve- 
locity disp ersion and maxi mum circular velocity are re- 
lated as in iFerraresd J2002) so that 



log(M BH ) = 5.051og(V 

1! 



4.41 



(4) 



which should be valid over the range Ub ~ 70 — 350km 
s _1 . We allowed different values of braking (77), rang- 
ing from none (77 = 1) to strong (77 = 0.1) which are 
shown in Fig. [lib). Clearly the braking parameter de- 
termines when the average mass captured reaches the 
total cloud mass. Our model predicts that the Edding- 
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ton ratio will be greatest near a critical value of V max , 
where the average mass captured reaches the total cloud 
mass. Evidently a nuclear strike on a black hole with 
mass < 1O 7 M0 yields much less mass available for accre- 
tion onto the black hole, even for large braking (r) = 0.1). 
However, a much larger fraction of M c i could be captured 
if the infalling cloud fragments into a stringy cloud com- 
plex in the galactic halo, effectively reducing r c i for a 
given Vmax- Cloud material not captured by the black 
hole will instead be gravitationally captured by the bulge. 
In this case, the material will either persist as part of the 
bulge ISM or will go down the route of star formation 
(with some unknown level of mixing) . This is consistent 
with observations of nuclear activity due to H\\ regions 
predominating in galaxies with smalle r bulges a nd there- 
fore smaller mass central black holes (|Hdl2008T ). 

6. ACTIVITY IN GALACTIC NUCLEI 

What sort of nuclear activity will result from a WHC 
impact? Once moderately dense, low angular momentum 
fragments of the WHC arrive in the galactic nucleus, star 
formation is likely, although dependent on poorly con- 
strained variables such as cooling rate and mixing with 
nuclear gas. The arrival of a large quantity of infalling 
gas may also initiat e activity via shocks or gra dual ac- 
celeration of GMCs ([Hopkins fc Hern auist 2003). In our 
model, shocked nuclear GMCs or the highest density, 
shocked fragments resulting from a WHC impact could 
generat e a nuclear OB ass ociation within ~ lOMyr after 
impact ([Clark et al.ll2005[) . A direct hit by a WHC frag- 
ment on the central black hole may result in a phase of 
lo w luminosity Bondi accr e tion. 

Kauffmann & Hcckman (2009) find two distributions 
of Eddington ratios among active galactic nuclei in the 
local Universe. Their first, a lognormal distribution, 
centered around a few percent of Eddington, is associ- 
ated with nuclear star formation, dominates among black 
holes < 10 8 M Q and depends critically on feedback from 
luminous accretion. The second distribution is a power- 
law which occurs in galaxy bulges where there is little or 
no ongoing star formation. We will only model the sec- 
ond mode since the first mode has a complex association 
with star formation, stellar evolution and AGN feedback. 
Although the distribution of dM/dt in this second accre- 
tion mode is a powerlaw, we will use the average value 
from our ACDM calculations to get an average activity 
fraction. We calculate the average fraction of galaxies 
of a given V max that are active for impactors of mass 
10 4 ~ 6 Mq and a range of braking and cloud densities. 
The lifetime of Bondi accretion is tb — M cap t/ (dM/dt) 
where 

dM/dt a lO- 3 nMf(^ 2 oo)- 3 M /yr. (5) 

Here n is the average density of the accreting gas, M 8 
is the black hole mass in units of 10 8 M Q and V200 is an 
average combination of gas spe ed and gas sound speed 
in units of 200 km/s (|Hol 12008). The average gas den- 
sity around the black hole is fairly unconstrained in our 
model and clearly there is trade off between more massive 
clouds that take longer to consume and higher gas den- 
sities which are consumed more quickly. Nevertheless, in 
this model the lowest average dM/dt is ~ lO _5 M0/yr 
around a 10 8 Mq black hole, assuming rj, V200 ~ 1 and 
n ~ 0.01cm -3 . In Fig. [ljc) we show some average 



values of braking and gas density that yield reasonable 
models where on average 1/3 of galactic nuclei are ac- 
tive. For 10 4 Mq clouds this requires n = 0.004cm -3 , 
while for 10 6 Mq clouds with some braking this requires 
n = 0.03cm~ 3 , and with no braking, n = 0.01cm -3 . 
Re call we are modelling t h e pow erlaw accretion mode 
of iKauffmann fc Heckmanl (j2009f ). which is significant 
among black holes of mass > 10 8 M©. Therefore, from 
Fig.[Hb),(c), our model requires 77 > 0.5, 0.3 for impactor 
clouds of masses 10 4 M Q and 10 6 M Q respectively. Note 
that in all cases the plot shows a characteristic shape 
where the average activity first increases as the amount 
of the cloud captured increases and then decreases as 
the accretion rate increases for the same amount of cap- 
tured cloud. Our model predicts that radiatively ineffi- 
cient Bondi accretion should be more commonly detected 
around larger mass black holes, which is consistent with 
observations of LINER activity in the local Universe (e.g. 
i Schawinski et al.ll2007t lHoll2008t IKauffmann fe Heckmanl 
2009). Future observations of low luminosity activity 
in galactic nuclei as a function of V max will be able to 
rule out certain parameters. Of course if r\ or n de- 
pend on V m ax then the shapes of the curves in Fig. QJc) 
might differ. MHD simulations of cloud fragmentation 
in the nucleus are required to establish whether the dis- 
tribution of densities and braki ng can in fact generate a 
powe rlaw distribution of dM/dt ((Kauffma nn fc Heckmanl 
2009) and we intend to do this in future work. Note 
that our model of Bondi accretion is unlikely to apply 
for M capt < M c i since the dominant activity in that 
case would be nuclear star formation, which is consis- 
tent with ob servations of smalle r mass central black holes 
(|Hol 120081 ). iCiotti fc Ostrikerl (|2007f l argue that lumi- 
nous accretion can occur in the aftermath of a nuclear 
starburst. This may yield the lognormal distribution of 
Eddington ratio observed in a fraction of galactic nuclei 
(jKauffmann fc Heckmanl l2009T l . but establishing this is 
far beyond the scope of the present work. 

We expect the number of WHCs to increase with red- 
shift as the gas accretion rate onto galaxies increases. We 
naively expect the number of WHC clouds and therefore 
the average cloud impact rate to scale with star formation 
rate or cx (1 + z) 31 from z — to z = 1 (Do hertv et al.l 
2006) , corresponding to a rate of ~ 1/ 30Myr for Milky- 
Way sized galaxies at z ~ 1. On these timescales, lumi- 
nous nuclear activity could result if radiatively efficient 
accretion occurs onto the central supermassive black 
hole. A stron gly decreasing Eddin gton ratio ove r cos- 
mic time (e .g. IVittorini et al J 120051 : IShankar et~aTl 120091 : 
iShenl l2009f ) could naturally produce a high duty cycle 
at z ~ 0, consistent with Fig. QIc). At even earlier 
epochs (z > 2.5), the bulk of the growth of supermassive 
black holes occurs via mergers in massive protoga lactic 
halos (~ 10 U - 13 M Q ) (ICavaliere fc Vittorinil[2002l) . At 
that time, gas accretion is primarily on random trajec- 
tories, so nuclear fueling by WHC impact would be at 
its maximum. A natural bridge can th en be made to the 
time before hot h alos are established ([Keres "et~alll2005l: 
iDekel et al.l [20091) . when gas accretion is all in the form 
of infalling gas. In the present epoch, our model predicts 
that clustering of activity should be random for similar 
sized galaxies (comparable V max ), since cloud impacts 
are random. This is consistent with observations of clus- 
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Fig. 1. — Plots of: (a) Average WHC impact rate (number/Gyr) 
on galactic nuclei versus maximum halo circular velocity (km/s) 
for impactor cloud masses of 10 4 (red),10 5 (green), 1O 6 M (blue). 
The corresponding black hole mass is shown on the top axis, (b) 
Average mass of 1O 4 M0 (red) and 10 6 Af o( j ot (blue) impacting 
clouds captured by the central black hole versus maximum halo 
circular velocity (km/s) for a range of braking (r) = 0.1, dashed 
lines;?} = 0.3, dot-dash lines, r\ = 1.0, solid lines), (c) Average 
fraction of activity in nuclei due to radiatively inefficient accretion 
versus maximum halo circular velocity (km/s) for 10 4 Mq (red) and 
10 6 Mq (blue) impacting clouds for a range of braking (r\ = 0.1 — 1) 
and a range of cloud densities (see text for details) . T he horizon tal 
dotted line corresponds to an activity fraction of 1/3 (Ho 2008). 



tering in LINERs IConstantin fc Vogelevl (|2006ft . 

7. CONCLUSIONS 

In this letter we introduce a new mechanism that deliv- 
ers large quantities of gas to galactic nuclei on astrophys- 
ically interesting timescales. We show that a single warm 
halo cloud (WHC) impact on a galactic bulge could po- 
tentially deliver a large mass (~ 1O 4_6 M0) of gas to the 
central regions of a galactic nucleus, in a singular event. 
Although there are considerable uncertainties in the pa- 
rameters of warm halo clouds, some representative num- 
bers suggest that the impacts occur on astrophysically 
interesting timescales and at a rate that must account 
for some or even all of the low luminosity activity ob- 
served in galactic nuclei in the local Universe. Based on 
analytic ACDM calculations of cooling halos, our model 
predicts an impact rate relatively independent of galaxy 
mass. Our model also predicts that larger mass black 
holes (> 10 7 M©) can capture significant fractions of the 
impacting cloud mass at some critical value depending 
on the nuclear braking (77). At this critical value, our 
model predicts the highest Eddington ratios. Below this 
critical mass, our model predicts that most of the cloud 
mass will not accrete onto the black hole. Instead, this 
material will mix with the ISM in the nuclear bulge or in- 
duce star formation, which may lead to delayed episodes 
of high Eddington ratio accretion. Finally, our model 
predicts that, for a reasonable range of cloud masses, 
densities and braking, the fraction of supermassive black 
holes accreting at very low Eddington ratios in the lo- 
cal Universe could be around ~ 1/3, which agrees with 
observations (|Hdl2008| ). 
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